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We generalize the de Almeida-Thouless line for the many-body Ising spin glass to the 
microcanonical ensemble and show that it coincides with the canonical one. This enables 
\^ • us to draw a complete microcanonical phase diagram of this model. 
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1. Introduction 

Spin glasses are magnetic systems in which spin interactions take random values 

'^ ! and remain fixed while the value of the spins may fluctuate. This is called quenched 

O ■ randomness. In such systems there occurs under certain conditions a thermodynamic 

phase transition where spins are frozen in random orientations, and thus their overall 

^ ■ configuration does not change, the spin-glass phase. 

^ i To calculate physical observables under quenched randomness we average over the 

random interactions. This is possible since extensive quantities, like the free energy, are 
t^^ [ self-averaging which means that they are the same for any realization of the random vari- 

ables. This usually entails averaging over complicated expressions, since the free energy 

for example, is the logarithm of the partition function. To deal with these difficulties, 
• 5— I ■ 
^ ■ the replica trick was invented. It relies on the simple equation: Inx = lim(x" — l)/n. 

'\J ■ n-i-O 

c^ I The notion behind this formalism is that it is much easier to evaluate some power of 

the partition function than its logarithm. 

In the theory of spin glasses, the de Almeida-Thouless (AT)-line^ is a boundary be- 
low which the replica symmetric solution of the Sherrington-Kirkpatrick (SK) model^ 
becomes unstable. It was generalized later by Gardner to include many-body spin glasses 
and replica-symmetry breaking.^ We established in a recent paper that for many-body 
spin glasses, the boundary between the paramagnetic and ferromagnetic phases is de- 
pendent on whether it was derived in the canonical or microcanonical ensemble.^ There 
is, however, no ensemble inequivalence where the spin-glass phase is concerned. It is 
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therefore expected that the result for the AT hne there is the same in both ensembles. 
In this paper we prove this assumption. 

Ensemble inequivalence occurs in systems with long-range interactions which are not 
additive. That means that when two subsystems with energy E are brought into contact 
their total energy will not be 2E in general. Many, at first glance, counter-intuitive 
effects appear in the microcanonical treatment of long-range interacting systems, like 
the appearance of negative specific heat. For a review of ensemble inequivalence, see 
e.g. Campa et al..^ 

The second purpose of this article is to draw the complete microcanonical phase dia- 
gram, including the AT line, minimal attainable energy, and the Nishimori line (NL)''''® 
of the many-body Ising spin glass. 

This paper is structured as follows. Following this introduction, in sec. 2 we introduce 
briefly the model. In sec. 3, we derive the microcanonical AT line, and discuss the 
condition of the NL in sec. 4. In sec. 5, we draw the complete phase diagram of the 
three-body spin glass. Paragraph 6 is devoted to concluding remarks. 

2. Model 

We restrict ourselves to the analysis of the Ising model with p-hodj, infinite-range 
interactions, called the infinite-range model. It is given by the Hamiltonian 

H = - ^ Jn..ipSi^...Si^ with Si^=±l, (1) 

ii<..<ip 

where the bonds are random numbers drawn independently from the distribution 

with mean Jqp\/Np~^, and jo is called the ferromagnetic bias. 

Using the replica trick, the microcanonical entropy density s, times the number of 
replicas n, is easily obtained.^ It is written with respect to spin-glass order parameter 

qai3 ~ X^i ^i^i/^ 2ind magnetization m^ ~ ^^ S^'/N as 

ns = -'^eaiQ'^)af3<^i3 -^QapQap - ^ m„ma -^ lu Tr 6^ , 

a/3 a</3 a 

^ = $^ga;3^"5^ + 5^m«5" (3) 

where ec = e -|- jomP^, Qap = Sa/^ + q^^^p^ and qaa = 0. The parameters qa/s and rha are 
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obtained from the saddle point conditions ds/dqap = and ds/dnia = as 

qal3 = 2pq^^'p^ Yl ^aQaaQj^l3, (4) 

m„ = 2pmP^-^jo Y^ Ql^ptp. (5) 

/3 

3. Microcanonical AT Line 

In this section we derive the microcanonical AT line for the replica symmetric (RS) 
and first-step replica symmetry breaking (IRSB) cases. 

3.1 General Considerations 

The entropy as given in eq. 3 cannot be solved in full generality. Usually, one makes 
some allowance for the symmetry in the replicas. For example, in the replica symmetric 
ansatz, one assumes that all replicas are identical and share the same values of order 
parameters. As it turns out, this ansatz leads to negative entropy at some finite tem- 
perature and is thus unphysical. However, it is valid at higher temperatures. Similarly, 
replica-symmetry braking ansatz may yield negative entropies at low temperatures. 

Let us now assume that we imposed some symmetry on the replicas and label the 
appearing parameters and quantities with a zero, i.e. the entropy s — )■ sq or the spin- 
glass order parameter q^o. To investigate the stability of this ansatz, we expand the 
entropy around the imposed solution in fiuctuations {qap — )■ q^p+dqap, iria — )■ m^^+dma) 
to second order s — )■ sq + X] dqapAap-ysdq^s + . . .. Since we used the saddle-point method 
to obtain the entropy from the sum of states, we require that the Hessian A has only 
positive eigenvalues. If this is violated at some point, then the imposed symmetry is 
obviously no longer valid. 

Calculation of the eigenvalues of the Hessian yields that there is indeed an eigenvalue, 
called the replicon mode, which may become negative at some temperature. However, 
it does not involve fluctuations of the magnetization. A careful examination of the 
expansion of the microcanonical entropy allows us to draw parallels to the canonical 
case and conclude that the microcanonical replicon mode will also be independent of 
fiuctuations of m. Therefore, we set subsequently diria = 0, which also implies drha = 0. 

Using the matrix expansion of the inverse of Q 

Q~^ ~ (Qo + dQ + d^QY^ ^ Qq^ - Qo\dQ + d^Q)Q^^ + Q^^dQQ^^dQQo', (6) 
and the expansion of the Fourier mode of spin-glass order parameter g^/? -^ Qai3 + '^^a/3 + 
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(Pqa/3, where 

dqafs = d f 2pql~^ ^ eaQ~^Qje,3 J 

^ 2p{p - l)ql-/eH'dq^^ - 2pql~'eH' ^ {{Qo%. + (Qo"')/3p) ^^^9.^,(7) 

with t = '^aiQo^)oii3i the expansion of the entropy in eq. (3) can be calculated in a 
straightforward manner. It is given by 

ns ^ ^"So - 2 X^ dqapdqap 

a<l3 

+^5^5^rfga/3rf'?.p((5°5^5'^50 - {S^S^){S'^S'^)). (8) 

a<j3 cr<p 

The average over the spins is the weighted trace, (■ ■ ■ ) = Tr(- ■ ■ )e^'-\ where 

We have already dropped the terms linear in fluctuations since the entropy must be 
extremized with respect to all its variables and thus, those terms disappear. 

If we take the n{n — l)/2 distinct elements of dq^fi and arrange them as a vector 
dqi, we can express the entropy as 

^ n(n-l)/2 

ns = nsQ- - ^ dqiAijdqj, (10) 

ij 

where the Hessian is made up of two matrices A = T^^ — G. Here, T is defined as the 
relation between dq^js and dq^p as given in eq. (7), i.e. dqi = '^jTijdqj, and the matrix 
G is implicitly defined in eq. (8). 

The calculation of the eigenvalues of A is somewhat tedious and we refer to the 
literature^ for details. It is sufficient for our purposes to note that the replicon mode 
can be expressed as 

\3={P-2k + R)-' + {P-2K + R), (11) 



where 



P = 2p{p - l)eH'q^-' - 2p'eHY^-'2{{Qo)-l + (QoV^), (12) 

K = -2phHY^~'2{{Qo)~i + 3{Qo)-;), (13) 

R = -2phHV'-'2{QorJ, (14) 
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are elements of the matrix T and 

P = 1-{S''S'^){S''SP), (15) 

K = (S^'S^) - (S^'Sf^) {S'^S") , (16) 

R = {S^'S^^S^'SP) - {S''S^){S''SP), (17) 

are elements of G. 

3.2 Replica Symmetry 

In the RS ansatz it is assumed that the values of the order parameters are equal in 

all replicas, m^ = m and q^p = q, for any a, (3. In this formulation, we have Q^^ = q^ 

ioi a ^ (5 and the components of A3 are easily listed as 

P-2k^R = Mp-l)f-He + jomn 

{1-qP) ^ ' 

P-2K + R = 1 - g2 _ 2(g - g2) + (S^S^S'^S^) - q^ 

= I Ducosh~'^{\/lqu + ih), (19) 

where Du = {2n)^^^'^ exp(—u'^ /2)du and we have used t = YlaiQ'^)aS ~ V(-'- + {n — 
l)g^) — > 1/(1 — g**). The condition that A3 is positive is given explicitly by 

2p(p-iyf-C+Mn'r ~" /^'■^°^''-^v^. + A). (20) 

This is the stability condition for the replica symmetric solution in the microcanonical 
ensemble. 

Then, we compare the result with the canonical one. The inverse temperature is 
defined as /3 = ds/de and is given in the RS case as /3 = — 2(e + jomP)/{l — q^). 
Inserting this relation into eq. (20), we see that, as far as the replica symmetric solution 
is concerned, the microcanonical and canonical AT lines coincide formally. 

3.3 First-Step Replica Symmetry Breaking 

With only replica symmetry, the spin-glass phase does not show up for p > 2, and 
we need to consider replica symmetry breaking. The IRSB ansatz is characterized by 
rUa = rn and g^/j = qo + {(li—(lo)Vai3(yX)—qi6ai3, where rja/six) is unity around the diagonal 
in blocks of size x and zero else. In the absence of the external field, solution of the 
IRSB saddle-point equations yields that the stable solution is always characterized by 
go = 0.^ Subsequently, we will set go = and relabel gi = g. 

Proceeding along the same line as the RS case, we obtain the IRSB AT condition 
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as 

{1 — {1 — x)qPy J Du[cosh{y/qu + m)]'^^'^ 

2p{p - l)qP-^{e + jomPy ^ J Du[cosh{^/qu + m)Y 
Using the IRSB energy-temperature relation (3 = — 2(e + jom^)/{l — (1 — x)q^), we see 
that eq. (21) is equivalent to the canonical AT condition.^ 

4. Nishimori Line 

In this section, we show the ensemble equivalence of the NL. In the original deriva- 
tion of the NL in the microcanonical ensemble,^ the random average is denoted as 

s = C H y"rfJ(.,...,)Po(Jn-.J^f^'Bjo- Yl 4-v lnTr5(i?-i/), (22) 

{ii---ip) \ {h---ip) ) 

where -Po(<^) is given by eq. (2) with jo = and C is an irrelevant constant. This 
is different from the standard measure (2). The crucial difference is that the gauge 
invariance can be utilized for eq. (22) and not for eq. (2). If we use eq. (22), the 
microcanonical version of the gauge invariant condition is given by e = —jo- This 
corresponds to the condition for the canonical case /3 = 2jo/ which implies that the 
difference between their distributions is irrelevant. We show explicitly this argument. 

If we use eq. (22) for the distribution, we obtain the same form as eq. (3) with the 
replacement 

Qa/3 -^ Qo.fi - ■mFoc^^fi- (23) 

However, this change does not affect the final form of the entropy. If we impose the 
IRSB ansatz, we obtain the replacement 

This difference disappears when we consider the n — )■ limit. Therefore, we conclude 
that the difference between eqs. (2) and (22) is irrelevant and we can utilize the gauge 
invariance on the NL e = —jo in our model. 

5. Phase Diagram 

Having obtained the microcanonical AT line in the previous section, we are finally 
able to draw the complete microcanonical phase diagram. In fig. 1 we compare the a) 
canonical and b) microcanonical phase diagrams. The canonical case was first obtained 
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Fig. 1. (Color online) Phase diagram of the many-body Ising spin glass with p=3 in the a) canonical 
and b) microcanonical ensemble. The thermodynamic phase boundaries are drawn in solid black lines, 
while the dynamical transition is drawn in solid red. The limit of the metastability of the ferromagnetic 
phase (spinodal line) is drawn in black dashed, while the AT line, below which the replica-symmetric 
solution is unstable, is drawn black dash-dotted. The blue dashed line is the NL. The limit of the 
stability of the IRSB solution of the spin-glass phase is drawn black dotted. 



in^° and the microcanonical phase diagram, except the AT hne, was drawn in.^ We 
re-draw here both for the sake of completeness. For jo < Jo^(~ 0.767) there is, in 
both ensembles, a horizontal, second-order phase boundary between a paramagnetic (P) 
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Fig. 2. (Color online) Microcanonical phase diagram in the (jo, e)-plane. The boundary between 
paramagnetic and ferromagnetic phases is drawn in black with circles. The spinodal line is drawn 
black dashed, while the AT line is shown black dash-dotted. The solid red line marks the dynamical 
transition. The full-RSB spin-glass phase exists only in the very narrow shaded region. The green 
dash-dotted line marks the minimal attainable energy. The blue dashed line is the NL. 
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and a IRSB spin-glass (SG) phase at Tc ~ 0.654. The stabihty-boundary is ensemble- 
equivalent, i.e. the SG phase becomes unstable at x ~ 0.335, Tb ~ 0.245 (black dotted), 
in the microcanonical as well as canonical ensembles. However, before the equilibrium 
P-SG transition takes place, there is a dynamical transition at T = Td(~ 0.686) (red), 
where the free energy develops an exponential number of minimal and the ergodicity 
breaks. ^^ The dynamical transition occurs at the same point in both ensembles. For 
Jo > Jo J there exists a replica-symmetric ferromagnetic (F) phase, which is separated 
from the P phase by a first-order transition. In the canonical ensemble this is a simple 
line (black), while in the microcanonical ensemble there is a region of phase coexistence, 
P-l-F in between.^ The ferromagnetic meta-stable states extend until Jq ~ 0.731, and 
the spinodal lines are shown black dashed. In both ensembles the ferromagnetic, RS 
solutions become unstable below the AT line, shown as black dash-dotted. Below this 
line there is a mixed (M) phase, where there is ferromagnetic order as well as RSB. The 
AT line is the same in both ensembles. The NL, with T = 1/2 Jq, is shown blue dashed. 
In fig. 2 we show the microcanonical phase diagram in the (jo, e)-plane. The F and P 
phases are separated, for jo > Jo') by a single line, shown black with circles. The replica 
symmetric F phase becomes unstable below the AT line, shown in black dash-dotted. 
The dynamical transition is at ed ~ —0.732 drawn in red. The condition for the NL 
reads e = —jo and is shown in blue dashed. The SG phase is stable for jo < Jq between 
ec(~ —0.769) and eb(~ —0.809). To estimate the value of the minimal attainable energy, 
emin, shown green dash-dotted, of the system we can take the energy value where the 
IRSB solution freezes. The energy where the entropy of the full RSB solution becomes 
zero at T = 0, lies at a higher energy, due to a peculiarity of the replica trick. Namely, 
the requirement that entropy be minimal with respect to the spin-glass order parameter 
and the RSB-parameter. A suggestive reason for this occurrence is that the term in the 
entropy as given in eq. (3) 

hm - V" qapQafi — > lim - I n^qoqo + -x^{qiqi - qoqo) - nqiqi \ (25) 

a</3 

changes sign at n = 1 since qo < qi- 

In fig. 3 we show the entropy as a function of the energy for fixed values of x 
at jo = 0.5. The equilibrium value of s is determined by taking the lowest available 
entropy value. Comparing the curves at constant x, we see that the energy where the 
entropy is zero increases when x decreases. For IRSB, however, this procedure ends 
at a; ^ 0.2. Lower values of x do minimize the entropy but lead to negative values. 
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Fig. 3. (Color online) Microcanonical IRSB entropy at jo = 0.5 for various values of x. At the point 
where the entropy is zero, indicated by an arrow, its derivative is less than infinity. This results in a 
temperature larger than zero at this point. 



We see furthermore that the energy where the entropy and temperature are both zero 
simultaneously when considering full RSB must lie at higher energy than emin- A similar 
statement holds for jo > Jo , for the ferromagnetic phase: The minimal attainable energy 
of the full RSB solution is expected to lie above emin which is the minimal energy of the 
IRSB solution. 

Between emin and the AT line there exists again, for jo > Jo^, a phase M which shows 
both, ferromagnetic order and RSB. On the other hand, the maximal attainable energy 
of the system is at e = 0. 

6. Conclusion 

In conclusion, we have shown that for Ising spin glasses with many-body interactions 
the AT condition yields the same curve in the canonical and microcanonical ensembles. 
Our significant result show that there is no ensemble inequivalence on the AT line. 
Since the ensembles are equivalent on NL as well and the AT line lies strictly at lower 
temperatures than the NL, we can surmise that there is no ensemble inequivalence below 
the AT line. This hypothesis is supported by another recent result^^ for spin glasses with 
integer spins, where the spin-glass phase transition can be ensemble inequivalent. There 
the AT line, however, terminates well before there is ensemble inequivalence. 
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